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ABSTRACT 


In  order  to  obtain  the  dynamical  structure  of  the  tropical 
atmosphere  in  equatorial  latitudes  a  short  range  prediction  experiment 
is  proposed.  The  initial  field  contains  the  intertropical  convergence 
zone  and  associated  disturbances  over  western  Pacific  Ocean  during 
March  1965*  A  complete  initial  state  is  constructed  using  a  consistent 
balance  system  of  equations .  The  procedure  involves  construction  of 
pressure,  temperature  and  vertical  motion  distributions  starting  from 
an  observed  rotational  part  of  the  wind  field.  It  is  shown  that  this 
procedure,  whose  validity  assumes  a  small  Rossby  number,  does  not  yield 
a  realistic  structure  of  the  vertical  motion.  A  short  range  prediction 
with  a  multi-level  prediction  model  yields  some  interesting  solutions 
in  the  vicinity  of  the  intertropical  convergence  zone.  During  the 
first  18  hours  of  prediction  an  adjustment  of  the  motion  and  the  mass 
field  ensues  with  gravity  inertia  oscillations.  A  detailed  discussion 
of  some  dynamical  aspects  at  2b  hours  is  presented.  An  important 
feature  of  the  model  is  a  parameterization  of  cumulus  convection  as  a 
function  of  large-scale  moisture  convergence.  The  role  of  cumulus 
scale  heating  in  the  vicinity  of  the  intertropical  convergence  zone 
is  explored  by  carrying  out  experiments  with  and  without  diabatic 


heating. 


List  of  symbols 


a 

fraction  (or  percent)  of  synoptic  scale  area  covered  by 

convective  clouds 

CD 

drag  coefficient 

C 

P 

specific  heat  of  air  at  constant  pressure 

e 

vapor  pressure  of  air.  Subscript  to  refers  to  water  surface 

Ev 

evaporation  from  sea  surface 

f’  fo 

Coriolis  parameter,  and  a  mean  value  of  the  Coriolis  parameter 

( respectively) 

F  ,  F 
x  y 

friction  force  per  unit  mass  of  air 

Fs 

flux  of  sensible  heat  from  ocean  to  the  atmosphere 

CD 

& 

Diffusion  of  heat  and  moisture,  respectively 

g 

acceleration  of  gravity 

H,  M 

heat  and  moisture  sources  and  sinks,  respectively 

V  hl 

sensible  and  latent  heat  (respectively)  per  unit  mass  of  air 

h 

terrain  height  above  sea  level 

,  j|,  Ik  unit  vectors  along  zonal,  meridional  and  vertical  directions 


i 

net  moisture  convergence  in  a  unit  vertical  column  extending 

from  900  millibars  to  the  top  of  the  atmosphere  (100  mb) 

ITCZ 

intertropical  convergence  zone 

eddy  kinetic  energy 

Kz 

kinetic  energy  of  mean  zonal  flow 

l 

number  of  grid  points  along  x-axis 

L 

latent  heat  of  vaporization  of  air 

m 

map  projection  factor 

P,  PQ 

pressure  and  a  standard  pressure  (=  1000  mb),  respectively 

PE 

eddy  potential  energy 

eddy  potential  energy 


Q. 

Q, 


R 


q 

Ro 

t,  e 

T 

a 


TS’ 


U,  V,  0) 

V,  \V 


0) 


C 


0) 

T) 

mn 

x,y,p>t 


potential  energy  of  mean  zonal  flow 
any  arbitrary  variable 

amount  of  moisture  required  to  replace  an  ambient  moisture 
sounding  by  a  local  moist  adiabat 

amount  of  moisture  required  for  condensation  warming  to 

replace  an  ambient  temperature  sounding  by  a  local  moist  adiabat 

gas  constant 

specific  humidity 

Rossby  number 

temperature  and  potential  temperature  of  air,  respectively 
temperature  of  air  at  anemometer  level,  assumed  same  as  ^qqq 
in  this  study 

temperature  and  specific  humidity  (respectively)  of  a  parcel 
raised  vertically  above  the  900  millibar  surface  with  no 
lateral  mixing 
temperature  of  water 

zonal,  meridional  and  vertical  (p)  velocity 

total  horizontal  wind  speed  and  wind  vector,  respectively 

average  vertical  velocity  over  region  occupied  by  convective 

clouds  inside  a  synoptic-scale  elementary  area 

average  vertical  velocity  over  a  region  not  occupied  by 

convective  clouds  inside  a  synoptic-scale  elementary  area 

synoptic  scale  vertical  velocity 

coefficient  of  a  9-point  local  surface 

independent  variables 


Z 


height  of  constant  pressure  surface 


0 

e 

equivalent  potential  temperature 

U  ,  < 

horizontal  and  vertical  eddy  diffusion  coefficients 

T  ,  T 

X  y 

surface  frictional  stresses  along  the  x  and  y  axes 

Ap 

vertical  grid  distance 

geopotential  height  of  pressure  surface 

zonal  mean  geopotential  on  a  boundary  wall 

$ 

latitude 

(i,  J) 

a  reference  grid  point  on  the  mercator  map 

P 

density  of  air 

*  *1000 

the  subscript  1000  refers  to  the  value  at  the  1000  millibar 

surface 

<l»,  c 

streamfunction  and  relative  vorticity 

% 

non-divergent  part  of  the  wind 

0 

static  stability  parameter,  =  -  — 

absolute  vorticity 

6 

the  beta  parameter,  - 

n 

RT 

Exner  pressure,  =  — 

X 

divergent  part  of  the  wind  field,  a  velocity  potential 

(  " ) 

bar  refers  to  a  zonal  average 

(  )' 

primes  at  a  symbol  refer  to  departures  from  a  zonal  average 

1.  Introduction, 


In  the  middle  latitudes  the  formation  of  an  extratropical  storm 
is  generally  regarded  as  a  baroclinic  instability  problem.  The  primary 
energy  source  that  leads  to  large  kinetic  energy  of  extratropical  vor¬ 
tices  is  the  available  potential  energy.  Phillips  (1956)  discussed  the 
energetics  of  this  class  of  disturbances  in  his  numerical  general  circu¬ 
lation  experiment. 

Lack  of  understanding  of  the  types  of  instabilities  in  the  tropics 
has,  so  far,  prevented  the  use  of  well  defined  numerical  experiments  to 
simulate  the  development  of  disturbances  that  are  frequently  portrayed 
from  observations  on  weather  maps. 

Some  success  in  the  understanding  of  the  formation  and  maintenance 
of  a  tropical  storm  (or  hurricane)  has  been  achieved  through  the  works  of 
Ooyama  (1963),  Charney  and  Eliassen  (196*0,  Kuo  (1965)  and  several  others. 
The  primary  reason  for  the  success  of  this  work  has  been  due  to  a  recogni¬ 
tion  of  what  is  called  a  conditional  instability  of  the  unsaturated 
tropical  air  mass.  Unlike  the  baroclinic  instability  problem  where  dis¬ 
turbance  with  horizontal  wavelengths  of  the  order  of  1000  or  more  kilo¬ 
meters  have  a  maximum  growth  rate,  for  the  tropical  conditional  instabil¬ 
ity  the  maximum  growth  rate  is  found  for  horizontal  wavelengths  of  5  to 
10  kilometers  in  scale  (Kuo,  1961).  This  latter  scale  is  now  recognized 
as  the  cumulus  scale.  Most  active  tropical  disturbances  contain  this 
cumulus  scale  in  the  form  of  organized  bands  of  convection.  There  is  a 
certain  cooperation  of  the  large  and  small  scale  the  nature  of  which  is 
a  central  problem  in  tropical  meteorology.  The  entire  belt  of  tropics  is 
frequently  conditionally  unstable,  namely 
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(0e  is  the  equivalent  potential  temperature  and  p  is  pressure.)  This  is 
observed  from  an  examination  of  atmospheric  soundings  below  the  500  mb 
surface.  However,  convective  scale  as  manifested  by  satellite  and  air¬ 
craft  observations  of  clouds,  is  present  only  over  portions  of  the  tropics. 
Thus,  there  is  evidently  a  supression  of  the  convective  motion  over  large 
areas  of  the  tropics.  This  might  be  controlled  by  the  large  scale  atmos¬ 
pheric  dynamics. 

A  mere  examination  of  the  streamline  features  on  a  weather  map  does 
not,  in  general,  give  indications  of  where  one  should  expect  convective 
motion  to  be  present  or  absent.  Several  tropical  meteorologists  present 
large  quantities  of  weather  maps  showing  merely  the  field  of  wind  direction. 
The  atmosphere  being  considerably  more  complicated,  not  many  questions  can 
be  thus  either  asked  or  answered. 

Because  of  large  data  voids  in  the  tropics  the  age  of  discovery  of 
new  phenomena  is  not  yet  over.  As  an  example  of  this  one  may  cite  the 
discovery  of  the  climatological  mean  troughs  of  the  upper  tropical  tropo¬ 
sphere  over  the  Atlantic  and  Pacific  Oceans.  These  were  deduced  from  com¬ 
mercial  aircraft  reports  of  upper  winds  during  recent  observational 
studies  by  Dean  (1956)  and  Sadler  (1967).  These  papers  do  not  deal  with 
either  the  three-dimensional  structure  or  the  dynamics  of  these  climato¬ 
logical  features.  There  is  thus  no  wide  recognition  of  the  importance  of 
these  studies.  The  data  is  sufficient  to  portray  cyclonic  or  anticyclonic 
streamline  features  on  a  constant  pressure  surface,  but  it  is  generally 
insufficient  for  study  of  the  dynamical  structure.  Several  climatological 
features  of  this  kind  must  exist  over  the  remaining  vast  uncharted  oceanic 
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regions  of  the  tropics.  Through  use  of  satellite  cloud  photographs  it 
might  be  possible  to  locate  regions  where  such  disturbances  are  present. 

This  paper  presents  some  features  of  the  tropical  circulations 
for  one  day,  March  1,  19^5,  as  part  of  an  investigation  which  is  being 
continued  over  a  longer  period  of  time.  The  study  deals  with  the  dy¬ 
namics  of  the  intertropical  convergence  zone.  The  approach  is  to  con¬ 
struct  a  consistent  three-dimensional  distribution  of  motion,  mass, 
temperature  and  moisture  fields  through  use  of  diagnostic  and  prognostic 
numerical  models.  The  large-scale  horizontal  wind  speed  in  the  tropics 
is  of  the  order  of  10  meters/sec  while  vertical  motion,  w  ,  is  of  the 
order  l  cm/sec  (Baumhefner,  1968).  The  latter  may  be  expressed  by  a 
relation  of  the  form 

_  . 

a)  =  a  a)  +  (l-a)  a)  (l) 

c 

where  ,al  is  fraction  of  the  synoptic  scale  area  covered  by  convective 
clouds,  and  u)  are  respectively  vertical  motions  in  the  cloud  and  out¬ 
side  the  cloud  in  a  unit  area.  Theories  of  parameterization  of  the  small 
scale  convective  motion  as  a  function  of  a  larger  scale  motions  enable 
us  to  determine  the  parameter  faT  subject  to  a  number  of  approximations, 
to  can  be  estimated  from  solution  of  balance  equations  and  adequation  or 
through  numerical  integration  of  primitive  equations. 

No  description  of  the  three-dimensional  structure  of  a  weather 
disturbance  can  be  complete  without  a  distribution  of  the  motion  field 
in  the  three  dimensions.  Any  investigation  of  energetics,  development 
and  instability  relies  heavily  on  the  distribution  of  vertical  motion. 

For  phenomenon  of  small  Rossby  number  (Ho  <  1),  a  measure  of  the 
vertical  motion  may  be  obtained  from  a  solution  of  the  quasi-geostrophic 
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to-equation  or  a  more  general  balance  a)-equation  (Krishnamurti ,  1968a). 

In  investigations  of  phenomenon  of  low  latitudes  (like  the  intertropical 
convergence  zone)  where  Ro  *,1,  the  validity  of  the  quasi -geostrophic  or 
balance  m-equation  is  questionable.  Primitive  equation  forecast  models 
contain  no  approximation  (except  the  quasi-static  approximation,  if  the 
vertical  equation  of  motion  is  a  hydrostatic  balance)  and  hence  a  short 
range  forecast  (using  a  quasi-geostrophic  or  balanced  initial  state) 
might  yield  a  reasonable  solution  of  the  three-dimensional  motion  field 
which  would  be  valid  for  large  Rossby  number.  Generally,  such  an  unbal¬ 
anced  initial  state  is  expected  to  produce  gravity  inertia  oscillations. 
These  are  damped  oscillations  and  there  is  an  ensuing  adjustment  of  the 
dependent  variables  into  a  large  Rossby  number  disturbance.  Some  map 
features  of  the  initial  state  are  lost  as  both  the  wind  and  the  pressure 
field  undergo  a  degree  of  adjustment.  The  process  of  adjustment  of  wind 
and  pressure  in  low  latitudes  in  the  presence  of  the  smaller  convective 
scale  is  a  very  complicated  and  important  theoretical  problem.  Simple 
formulations  of  this  problem  by  Rosenthal  (1965)  and  Matsuno  (1966a)  are 
important  contributions. 

We  shall  show  in  this  paper,  that  the  numerical  solution  of  the 
primitive  equations  (starting  from  solutions  of  a  diagnostic  balance 
model)  undergoes  an  adjustment,  the  theoretical  nature  of  which  is  not 
quite  clear,  but  the  process  is  certainly  of  interest.  It  should  be 
noted  that  the  adjustment  and  the  damping  of  the  gravity  inertial  oscilla¬ 
tion  is  brought  about  by  the  diffusion  of  momentum  and  potential  tempera¬ 
ture.  The  diffusion  coefficients  are  selected  somewhat  arbitrarily  in 
these  studies*,  the  choice  was  made  by  trial  and  error.  In  our  study  after 
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a  period  of  approximately  12  hours  these  oscillations  become  very  small 
in  amplitude.  We  are  thus  able  to  present  a  description  of  the  three- 
dimensional  distribution  of  mass,  moisture  and  motion  field  for  large 
Rossby  number  phenomena  through  short  range  numerical  prediction  of 
atmospheric  variables.  Two  questions  that  come  up  in  the  light  of  these 
comments  regarding  the  oscillations  are: 

i)  How  good  is  the  initial  data  for  a  description  of  the  large 
scale  distribution  of  the  dependent  variables  in  regions  near  the  equa¬ 
tor  or  in  the  tropics? 

ii)  How  much  of  the  initial  map  features  is  lost  during  this 
process  of  adjustment? 

We  are  concerned  here  with  a  study  of  the  intertropical  converg¬ 
ence  zone.  There  is  a  considerable  amount  of  data  in  the  western  Pacific 
Ocean  to  describe  the  large-scale  flow  both  at  the  surface  level  and  up 
to  the  200  millibar  surface.  This  is  borne  out  by  a  number  of  descript¬ 
ive  studies  in  this  region  (Yanai,  196l).  This  data  is  supplemented  by 
a  number  of  aircraft  reports  (AIREPS)  that  enable  a  better  analysis  of 
the  700  and  250  millibar  charts.  The  initial  analysis  at  the  lowest 
levels  contains  features  like  a  convergence  line  (the  ITCZ).  At  the  upper 
level  features  like  warm  high  pressure  cells  are  found  to  the  north  and 
south  of  the  low  level  convergence  line.  The  anticyclonic  outflow  from 
these  upper  level  high  pressure  cells  (size  lOOOnds  of  kilometer  scale) 
is  an  easterly  wind  maximum  over  the  equatorial  latitudes  and  lies 
directly  above  the  surface  level  convergence  asymptotes.  The  damped 
oscillation  of  the  prediction  does  not  destroy  these  large-scale  features. 
Hence,  a  description  of  the  predicted  fields  may  be  considered  a  description 
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of  the  large-scale  dynamical  structure  of  the  ITC.  This  is  the  principal 
motivation  for  carrying  out  this  investigation. 

2.  A  primitive  Equation  Model 

(a)  Advective  scheme  and  basic  equations: 

The  main  features  of  the  advective  scheme  are  essentially  the 
same  as  those  proposed  by  the  author,  Krishnamurti ,  1962.  Similar  quasi- 
Lagrangian  schemes  have  been  proposed  by  Leith  ( 196U )  and  Okland  ( 19^5 ) - 
Leith  carried  out  long-term  numerical  integrations  in  a  general  circula¬ 
tions  model  while  Okiand  proposed  its  use  in  short-term  integrations  for 
a  two-level  model.  No  non-linear  advective  terms  (\V  •  VQ)  appear  explic¬ 
itly  in  this  formulation  and  thus  no  obvious  non-linear  computational 
instabilities  are  present.  In  a  recent  review  of  various  problems  gener¬ 
ally  encountered  in  numerical  weather  prediction,  Elsasser  (1968)  points 
out  that  (jkland’s  Lagrangian  advective  scheme  was  one  of  the  important 
factors  in  the  successful  prediction  of  cyclogenesis  over  short  time 
periods  of  the  order  of  two  days. 

The  main  outline  of  advective  scheme  will  be  discussed  here. 
Consider  the  basic  system  of  the  primitive  equations  expressed  in  the 
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This  may  be  regarded  as  a  closed  system  for  the  seven  variables  u,  v,  0, 

3  z 

^1000’  “  3p  °r  alternatively  a  closed  system  for  u,  v,  0,  q,  z  and 
a),  with  appropriate  boundary  conditions.  is  a  quasi-Lagrangian  opera¬ 


tor  given  by 
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Here  m  is  a  map  factor,  in  our  case  for  a  mercator  projection  expressed 
by  the  relation 

m  =  sin  85  cosec  (90  -  $)  (10) 

where  $  is  the  latitude,  m  =  1  at  $  =  15°N  and  15°S5  at  the  equator 
m  =  0.996.  Equations  (2)  and  (3)  are  the  momentum  equations, 

(U)  is  a  form  of  the  thermodynamic  energy  equation.  Equation  (5)  is  a 
statement  for  conservations  of  water  vapor.  Equation  (6)  is  a  tendency 
equation  of  the  pressure  frame  where  the  role  of  orography,  h,  is  in¬ 
cluded.  (7)  and  (8)  are,  respectively,  the  continuity  and  hydrostatic 
equations.  Heat  sources  and  sinks,  H,  and  the  moisture  sources  and 
sinks,  M,  are  discussed  later  in  the  paper.  F^  and  F  ,  FQ ,  and  F^  repre¬ 
sent  frictional  and  diffusive  transfers  of  momentum,  heat  and  moisture, 


and  are  also  discussed  later. 
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The  principle  of  the  quasi-Lagrangian  advective  scheme  is  to 
locate  the  past  position  (P)  of  a  parcel  on  a  pressure  surface  during 
each  time  step,  the  parcel  arrives  at  a  grid  point  0  (I,  J)  at  the  end 
of  each  time  step,  see  figure  1.  During  the  course  of  movement  from 
P  to  Q  non-conservative  changes  take  place  which  are  incorporated  in 
the  advected  properties. 

The  iterative  procedure  utilizes  a  local  surface  over  9  points 
of  figure  1  defined  by  the  relation 


2  2 


l  l 


n  X®  Yn 
mn 


(11) 


The  nine  coefficients  ri^  are  predetermined  for  each  point  and  do  not 
vary  during  the  course  of  the  integration.  X  and  Y  are  distances  in 
the  zonal  and  meridional  directions  between  the  central  point  I,  J  and 
the  other  8  points  of  figure  1. 

Distance  (PO.)  =  tyx2  +  Y2. 

Point  P  is  determined  by  a  predictor  corrector  successive  ap¬ 
proximation  procedure,  the  first  guess  is  defined  by  the  relations 

X  *  -  ' 


(12) 


Y  '  "  At 

The  variables  u,  v,  w,  z,  0,  q.  A,  B,  C,  D  and  E  are  determined  at  the 
point  P  using  relation  (ll)  and  the  values  of  X  and  Y  defined  in 
equation  (12). 

Dote  that  At  is  a  time  step:  the  subscript  t-At  is  the  previous  time 
period*  and  t  would  be  the  current  time  for  which  a  prediction  is 
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desired.  The 

first  guess  prediction  of  the 

dependent  variables  is 

carried  out  according  to  the  relations, 

u(Q)t 

= 

U<P)t-At  +  A<P)t-At  4t 

*\ 

v(Q)t 

= 

V<P>t~At  *  B(P)t-At  4t 

e(Q)t 

= 

e(P)t.4t  *  c<p)t_st  At 

\  (13) 

/ 

q(Q)t 

= 

■><P>t-At  *  D<P)t-At  4t 

Z1000(Q)t 

Z1000^P^t-At  +  E^P^t-At 

At 

) 

Z(«)t 

and 

gj(q)  are  determined  from 

u 

the 

diagnostic  relations 

(7)  and  (8). 

The  second  guess  defines  a  location 

of  the  point  P  through  the 

relations 

X  = 

- 

(“(P)t-At  *  4t/2 

N 

y  (im 

Y  = 

- 

fv(P)t-At  *  V<Q)t'  4t/2 

> 

Variables  u,  1 

7,  0) 

,  Z,  6,  and  q  are  interpolated 

at  point  P  (time  t-At) 

and  the  prediction  at  point  Q  is  carried  out  using  the  relations 

u(Q)t 

= 

u(P)t-At  +  A(Q)t  4t 

v(Q)t 

= 

-<P>t-At  *  B(9)t  4t 

e(o)t 

= 

0(P)t-At  *  C(Q>t  4t 

(15) 

q(Q)t 

- 

«<P>t-At  *  D(9)t  4t 

/ 

ziooo 

<o)t 

’  WP)t-At  *  E(«>t 

At 

10 


Z(Q)t  and  u)(Q)t  are  again  determined  from  equations  (7)  and  (8). 

This  two-step  predictor  corrector  is  similar  to  Matsuno’s  (1966b) 
time  differencing  scheme  for  Eulerian  and  spectral  equations  (  see  also 
Lilly  (1965)).  Terms  containing  vertical  advection,  pressure  gradients, 
Coriolis  force,  friction,  heating  and  other  non-conservative  features  on 
the  right  hand  side  of  equations  (2)  through  (7)  are  incorporated  in  A, 

B,  C,  D  and  E.  A  time  step  At  =  10  minutes  is  used  in  the  present  study. 

b)  Domain  of  integration,  distribution  of  variables: 

The  domain  of  integration  extends  from  10U°E  longitude  to  128°W 
longitude  over  the  Pacific.  Northern  and  southern  limits  of  the  domain 
extend  from  28°N  to  28°S.  The  region  between  152°W  to  128°W  is  used  to 
prescribe  a  cyclic  continuity  in  the  east-west  direction.  This  is 
carried  out  by  a  quadratic  interpolation  of  the  initial  analyzed  data 
such  that , 

Q(10U°E)  =  Q(128°U)  (16) 

where  Q  is  any  analyzed  variable.  It  must  be  noted  that  the  initial 
state  over  the  entire  domain  is  a  consistent  balanced  hydrostatic  set, 
cyclic  continuity  merely  removes  east-west  boundaries  in  both  the  dia¬ 
gnostic  and  prognostic  studies  presented  here. 

In  a  recent  comparative  study  of  several  finite  difference 
analogs  for  primitive  equations,  Grammeltvedt  (1968)  states  that  if  in¬ 
tegrations  are  performed  for  more  than  three  days  then  more  than  15  grid 
points  per  wavelength  are  needed  to  describe  with  accuracy  the  movement 
and  development  of  the  shortest  wave  which  initially  is  carrying  a  sig¬ 
nificant  part  of  the  energy.  This  is  an  important  result  for  numerical 
weather  prediction.  This  would  imply  that  for  studies  of  development  and 


11 


movement  of  a  disturbance  in  the  tropics,  whose  wavelength  is  of  the 
order  of  1000  km,  a  mesh  size  of  the  order  of  6 0  km  would  be  required. 

A  fine  mesh  model  will  be  widely  used,  when  faster  computers  with  larger 
core  storage  become  available.  Phase  speed  errors,  through  reduction  of 
truncation  errors,  can  be  reduced  appreciably  by  this  procedure.  In  the 
present  study  we  use  a  mesh  size  of  2°  latitude  by  2°  longitude  and  its 
limitations  are  thus  obvious.  In  certain  large-scale  general  circulation 
models,  mesh  sizes  of  the  order  of  5°  latitude  are  being  used. 


In  the  present  study  there  is  no  staggering  of  variables  in  the 


horizontal;  however,  there  is  a  staggering  of  dependent  variables  in  the 
vertical.  It  is  a  five-level  model  and  staggering  permits  all  nine  levels 
to  be  used.  This  vertical  structure  is  primarily  dictated  by  the  nature 
of  the  dynamical  and  thermodynamic  equations,  see  figure  2. 

c)  Boundary  conditions  for  the  primitive  equation  model 

As  mentioned  earlier,  there  are  no  boundaries  in  the  zonal  direc¬ 
tion  in  this  study.  Our  basic  grid  extends  from  28°N  to  28°S.  There  is, 
however,  an  implicit  north  and  south  pole  in  the  model.  At  the  north  and 
south  poles  we  assume  at  each  time  step 


ziooo  (North  pole)  =  l  Z1000  (28°N)/£ 

X 


6  (North  pole)  =  £  6  (28°N)/£, 


X 


q  (North  pole)  =  £  q  (28°N)/£ 


x 


(17) 


x 


0  (South  pole)  =  J  0  (28°S)/£ 


x 


q  (South  pole)  =  J  q  (28°S)/£ 


x 


12 


u,  v  and  to  at  the  two  poles  are  assumed  zero  during  the  integra¬ 
tion.  A  linear  interpolation  between  pole  and  28°  is  used  to  calculate 
the  values  at  30°N  and  30°S  during  each  time  step.  These  values  were 
used  as  boundary  conditions  during  the  course  of  the  integration.  This 
scheme  was  found  to  be  stable  and  values  of  the  vertical  motions,  the 
most  sensitive  parameter  of  this  model,  do  not  become  unreasonable.  This 
boundary  condition  does  not  in  any  way  represent  the  behavior  of  the  real 
atmosphere  at  30°M  and  30°S. 


d)  The  development  in  the  primitive  equation  model 

In  quasi-geostrophic  models  development  is  generally  discussed  in 
terms  of  processes  producing  individual  changes  of  vorticity  via  the  vor- 
ticity  equation.  These  changes  are  primarily  attributable  to  the  pres¬ 
ence  of  fields  of  divergence  and  convergence.  In  primitive  equation 
models  surface  development  may  be  described  via  the  follbwing  equations. 

rlOOO 

<ia> 

100 


“lOOO 


_D_ 

Dt 


(Z1000  -  h) 


“1000  pg  61000 


(18b) 


If  h  =  0  (no  terrain)  and  if  the  vertically  integrated  divergence 
is  zero,  then 


and 


“lOOO  0 


DZ1000 

Dt 


0 


(19) 


(20) 


i.e.,  parcels  moving  on  the  1000-millibar  surface  will  conserve  their 
heights  and  no  deepening  or  filling  of  surface  systems  is  possible. 
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Up  and  downslope  motions  over  large-scale  terrain  features  permit  de¬ 
velopment  of  surface  systems  such  as  upslope  anticyclogenesis  and  a 

lee  cyclogenesis.  Net  convergence  leads  to  a  sinking  motion  (a)^000>  ^ 

DZ 

and  — >  O)  >  while  the  converse  is  true  of  net  divergence.  One 
Dt 

would  think  that  since  the  vertically  integrated  divergence  is  a  small 
difference  between  large  numbers  (at  each  level)  the  results  of  such 
vertical  integration  will  be  very  sensitive  and  only  spurious  results 
can  be  obtained.  This,  however,  definitely  is  not  the  case.  A  bal¬ 
anced  initial  state  starts  out  with  minor  gravity  inertia  oscillations 
which  damp  rapidly  with  time  and  very  small  and  reasonable  magnitudes 
of  ^2000  are  produced  by  the  time  integrations,  see  figure  12. 

This  has  been  verified  by  nearly  all  of  the  research  workers  who  have 
integrated  the  primitive  equations.  A  synoptic  meteorologist  finds 
it  very  difficult  to  understand  this  feature  of  the  primitive  equation 
models  when  he  knows  that  he  has  only  a  limited  success  in  the  use  of 
the  so-called  tendency  equation.  The  apparent  reason  for  this  lies  in 
the  fact  that  a  primitive  equation  model  generates  a  dynamically  con¬ 
sistent  momentum  field  with  sufficient  accuracy  such  that  the  net  di¬ 
vergence  and  tendencies  are  small,  as  is  the  case  for  the  real  atmosphere. 
Direct  observations  do  not  yield  small  net  divergences  or  tendencies 
because  winds  usually  are  in  error.  A  detailed,  dynamic  explanation  of 
this  process  of  adjustment  would  be  of  considerable  synoptic  interest. 

In  the  vicinity  of  mountains  w^oOO  acquires  magnitudes  very  close 
to  that  given  by  the  lower  boundary  conditions  of  the  m-equation ,  i.e. 

^1000  =  ~  ^p1000  W1000  ‘  This,  again,  is  somewhat  surprising  because 

in  the  primitive  equation  model,  terrain  h  modifies  ^000  v^a  a  ra^her 
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complicated  cycle  of  operations. 

h  modifies  Z^qqq  via  equation  (6) 

^1000  ^e^erm:*-nes  Z  via  equation  (8) 

Z  determines  u,  v  via  equations  (2)  and  (3) 
u,  v  determines  ^000  v^a  equation  (7) 

The  reason  for  this  must  again  be  that  a  dynamically  consistent 
primitive  balance  is  indeed  obtained  when  the  initial  state  is  given  by 
the  balance  equations.  The  nature  of  such  a  primitive  balance  in  the 
vicinity  of  the  equator,  where  the  Rossby  number  is  greater  than  unity, 
has  not  been  explored  in  theory  as  yet. 

3.  Initialization 

A  primitive  equation  model  is  currently  being  used  in  operational 
weather  prediction  at  the  National  Meteorological  Center  in  the  United 
States  (Shuman  and  Hovermale,  1968).  The  initial  wind  field  is  obtained 
from  a  solution  of  the  so-called  linear  balance  equation.  This  relation¬ 
ship  between  the  wind  field  and  the  pressure  field  does  not  permit  any 
divergence  or  vertical  motions  at  the  initial  time.  A  similar  initial 
state  was  proposed  for  the  Bushby  and  Timpson's  (1967)  model  by  Benwell 
and  Bretherton  (1967)  in  a  recent  paper.  It  is  pointed  out  by  the  lat¬ 
ter  authors  that  this  procedure  suppresses  inertio-gravitational  oscilla¬ 
tions.  These  ideas  are  by  no  means  new,  Charney  (1955),  Hinklemann  (1951) 
and  Smagorinsky  et  al  (1965)  have  proposed  the  use  of  non-divergent  ini¬ 
tial  states  for  integration  of  primitive  equations.  Phillips  (i960),  on 
the  other  hand,  showed  that  an  initial  divergence  equal  to  that  prescribed 
by  quasi-geostrophic  dynamics  would  be  required  to  suppress  the  unwanted 
gravity  inertio-oscillations . 
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In  this  paper  we  propose  a  rather  detailed  diagnostic  balance 
model  to  define  a  consistent  initial  state.  Thus,  initially  both  diver¬ 
gence  and  vertical  motions  are  present.  A  part  of  our  interest  will  be 
to  describe  the  three-dimensional  synoptic  structure  and  the  dynamics  of 
the  motion,  mass,  thermal  and  moisture  variables  at  the  initial  state 
from  balance  equations.  This  study  departs  somewhat  from  earlier  study 
on  initialization  by  the  author  (Krishnamurti ,  1968a);  here  we  propose 
an  initialization  of  the  dependent  variables  across  the  equator  in  a 
tropical  belt.  This  poses  seme  problems  if  the  streamfunction  were  to 
be  obtained  from  observed  geopotential  heights,  but  it  is  a  relatively 
simple  matter  if  the  streamfunction  is  deduced  from  the  observed  wind 
field.  The  initialization  procedure  is  as  follows: 

(a)  Streamfunction  ^ 

The  nondivergent  part  of  the  wind  is  obtained  from  a  solution 
of  the  equations 

v2ip  =  c  (21) 

and  \V  =  |k  x  Vip  (22) 


where  £  is  the  relative  vorticity  of  the  observed  wind.  Boundary  condi¬ 
tions  for  i|>  are  obtained  by  the  procedure  outlined  by  Hawkins  and  Rosen¬ 
thal  (1965)  (version  2  of  their  work).  This  particular  choice  of  boundary 
condition  permits  no  net  mass  flux  from  the  domain;  however,  i|>  is  permitted 
to  vary  along  the  boundaries . 

(b)  Initial  distribution  of  geopotential  heights 

The  distribution  of  geopotential ,<J>  =  gZ,  is  obtained  by  solving 
the  balance  equation 


VI  = 


,  .  m  +  2J  (ft  ,  ft) 


(23) 
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Since  i()  (x,  y,  p)  is  known  the  right  hand  side  is  defined  in  the  domain 
and  the  solution  of  this  elliptic  equation  is  a  straightforward  problem. 
At  the  boundaries  we  assume  <j>  to  be  given  by  the  relation 

<{>B  =  4>  +  fQ  (\p  -  (23a) 

4>b  is  the  value  of  <j>  at  the  northern  or  southern  walls  at  28°N  or  28°S. 

$  is  the  zonal  mean  value  of  <p  obtained  from  observations  of  geopotential 

height  of  the  northern  or  southern  wall. 

-k  -1 

fg  =  10  sec  at  northern  wall, 

-if  -i 

and  fg  =  -  10  sec  at  southern  wall. 

\p  is  the  value  of  the  streamfunction  at  any  point  along  the  wall  and 
ip  the  zonal  mean  value  along  that  wall. 

The  initial  distribution  of  potential  temperature  is  given  by 


Potential  temperatures  computed  by  this  procedure  were  found  to  depart 
by  as  much  as  3  degrees  in  the  middle  latitudes  from  observations. 
Maximum  departures  were  of  the  order  of  5  degrees  at  lower  levels. 

This  difficulty  could  have  been  overcome  by  using  observed  tempera¬ 
tures  directly  in  the  model;  this,  however,  would  cause  the  sacrifice 
of  the  initial  balance  relationship  of  the  variables.  It  is  not 
a  priori  clear  how  serious  an  effect  this  would  have  had  on  the  initial 
amplitudes  of  gravity  inertio-oscillations . 

( c )  Divergent  part  of  the  wind,  X,  and  vertical  velocity 

A  detailed  analysis  of  the  system  of  equations  required  to  solve 
for  X  and  oj  is  presented  elsewhere  by  the  author  (Krishnamurti ,  1968a). 
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Here  we  shall  merely  present  the  equations  and  their  boundary  conditions. 
Omega  equation 


_2  ^  .2  32oj 

V  oto  +  f  — r- 

3p 


+  "’2  J<*-6)  +  k 


3r  3t 
3x  3  y 


-  2  —  ~  j rli  iih  _  f  _2_  fru^y) 

2  3t  3p  j^3x’  3y-*  f  3p  ^  X' 


r-p  ’2"  * f  ,7  ("  &  ^ 


*  r  F»  (7“  •  7  IS  -  f  £  f'x  •  vO 


ap 


3p 


(25) 


Continuity  equation 


v2x  = 

3p 


(26) 


Vorticity  equation 


V2  If  =  -  J(*.sa)  *  n 


V?  +  c  V  X 
a  a 


3  „2,  „  „  3if>  3 

-  o>  ^  V  *  -  Vm  •  V  -  -  g  - 


3F  3F 

_ y  x 

3x  3y 


(27) 


Boundary  conditions 


At  p  =  100  mb,  a)  =  0 


At  southern  wall,  y  =  y^,  to  =  X  =  |^-  =  0 


At  northern  wall ,  y  =  y  ,  uj  =  X  =  =  0 

C.  O  t 


At  p  =  1000  mb ,  to  =  -  g 


1000 


R  0 


1000 


j(i|),  h)  -  VX  .  Vh 


■] 
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The  solutions  for  w,  X  and  are  obtained  by  a  successive  approxi- 

<3t 

mation  procedure.  This  involves  three-dimensional  relaxation  for  m  and 
two-dimensional  relaxation  for  X  and  .  These  three  equations  are 
coupled.  Conditions  for  ellipticity  of  these  equations  are  discussed  by 
Krishnamurti  (1968a). 

and  are  frictional  forces  per  unit  mass  of  air,  along  the 
x  and  y  axis:,  H  is  net  diabatic  heating  per  unit  mass  of  air.  These  are 
discussed  in  some  detail  in  a  separate  section  of  this  paper.  The  same 
forms  of  friction  and  heating  are  retained  in  the  initialization  and  pre¬ 
diction  equations. 

The  total  initial  wind  is  given  by 


\V  = 


(28) 


Initial  distribution  of  moisture  is  a  scalar  hand  analysis  of 
relative  humidity  at  standard  meteorological  levels,  this  is  converted 
into  fields  of  specific  humidity  q  at  the  900,  TOO,  500  and  300  millibars 
surfaces.  The  closed  system  of  initial  dependent  variables  are  u,  v,  m, 
Z,  6  and  q,  figure  2  shows  the  vertical  staggering  of  variables. 


I4.  Diabatic  heating  and  parameterization  of  cumulus  scale 

In  the  present  model  diabatic  heating  per  unit  mass  of  air  is 
split  into  two  components. 

H  =  Hs  +  Hl  (29) 

where  H  is  the  transfer  of  sensible  heat  from  the  ocean: 

is  the  release  of  latent  heat. 

Ue  have  not  considered  radiative  effects. 
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(a)  Sensible  heat  flux: 

A  simple  flux  formula  (Jacobs,  1951),  was  used  to  define  this 
quantity. 


Fg  =  0.0329  V  (Tw  -  Tfl)  (30) 

The  unit  of  flux,  F  ,  is  mb  meter  sec”\  V  being  the  wind  speed 
in  meters/sec  at  anemometer  level  and  (t  -  T  )  the  temperature  differ¬ 
ence  between  water  at  the  ocean's  surface  and  air  at  the  anemometer 
level.  The  constant  of  proportionality  has  dimensions  mb/°C.  The  rate 
of  heating  Hg  may  be  expressed  by  the  relation 


H 


S 


(31) 


Potential  temperature  6  is  a  dependent  variable  and  it  appears 
at  the  900,  700,  500  and  300  millibar  surfaces.  It  is  assumed  that 
sensible  heat  transfer  modifies  the  potential  temperature  of  the  air  at 
the  900  millibar  surface,  and  the  flux,  Fg,  is  assumed  zero  at  800  mb 
and  above. 

We  express  sensible  heat  at  900  millibar  surface  by  the  formula 

HS  =  fefT  °-0329  V  <T„  -  Ta> 

=  0.00323  V  (T  -  T  )  (32) 

w  a 

2  -3 

Here  H  has  units  meter  sec  .  If  the  air-sea  temperature  difference 

is  of  the  order  of  1°C  and  the  wind  speed  is  of  the  order  10  mps,  then 

a  change  of  one  degree  in  potential  temperature  may  he  expected  over  a 

2U-hour  period.  This  would  require  a  representative  analysis  of  the 

initial  distribution  of  (T  -  T  )  over  the  oceans.  Because  of  general 

w  a 
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lack  of  data  only  a  rather  smooth  representation  of  this  field  is  pos¬ 
sible  at  this  time.  Radiative  cooling  is  a  serious  emission  of  the 
present  study  because  it  is  also  of  the  order  of  a  degree  centigrade 
per  day. 

(b)  Latent  heat 

Ht  may  again  be  divided  into  two  forms  of  latent  heat,  stable 
and  conditionally  unstable.  The  former  may  be  important  in  regions  of 
dynamic  ascent  of  moist  and  absolutely  stable  air.  The  latter  may  be 
important  if  there  are  regions  where  certain  prescribed  conditions  for 
the  evolution  of  cumulus  scale  convective  elements  are  met.  For  the 
present  we  have  omitted  the  stable  form  of  latent  heat  as  the  tropical 
air  mass  is  generally  known  to  be  conditionally  unstable  in  the  lower 
troposphere  and  our  interest  in  this  study  is  an  equatorial,  tropical 
belt . 


The  procedure  for  parameterization  of  the  cumulus  scale  as  a 
function  of  the  large-scale  motion  with  the  mutual  modification  of  the 
two  scales  is  carried  out  along  the  lines  of  Kuo's  (1965)  formulation 
for  a  tropical  storm. 

Heat  is  released  if  the  following  criteria  are  met: 

i)  _  —  <  0  at  any  level  above  1000  millibars.  This  measures 

3  P 

the  moist  conditional  instability. 

ii)  Net  moisture  convergence,  I  >  0,  in  vertical  columns.  The 
column  extends  from  the  top  of  the  friction  layer  to  the  top  of  atmos¬ 
phere  (100  mb). 


The  heating  function,  H  ,  is  defined  by  the  relation: 

L 


H 


L 


-  T) 


(33) 
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Where  (T  -  T)  is  the  temperature  difference  between  cloud  element  and 

D 


tropical  environment  r,  is  the  specific  heat  of  air  at  constant  pres¬ 
sure;  At  is  a  cloud  time  scale,  where  in  a  time  integration  problem  we 
use  At  =  the  time  step  of  the  primitive  equation  model  (Kuo,  1965);  and 
"a"  is  a  measure  of  the  fraction  of  the  grid  scale  area  covered  by  act¬ 
ive  convective  elements. 

The  calculation  of  "a"  involves  determining  the  following  three 
integrals . 

Net  moisture  convergence  above  friction  layer. 

W900  q900 


I  =  * 

G 


rPT 


q  V  dp  - 


900 
■PT 


g 


On 


gk  (a-s  -  «>  * 


(3U) 


(35) 


900 


*  C  ^  • 


(36) 


and 


a  =  I/(Q1  +  Q2) 


(37) 


T  and  q  are  respectively  the  temperature  and  specific  humidity 

O  D 

of  a  parcel  raised  dry  adiabatically  from  900  millibar  surface  to  the 
lifting  condensation  level  and  then  moist  adiabatically  to  a  level  p^ 
where  intersects  the  environmental  temperature.  A  detailed  analysis 
of  this  parameterization  procedure  and  aspects  of  calculation  are  dis¬ 
cussed  elsewhere  by  the  author  (Krishnamurti ,  1968b).  In  this  earlier 
study  diagnostic  calculations  of  vertical  motions  and  various  friction 
layer  parameters  were  evaluated.  It  was  shown  that  convective  banded 
motions  as  revealed  by  airborne  radar  and  that  calculation  of  vertically 
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integrated  moisture  convergence  and  the  field  of  the  parameter  "a" 
exhibit  reasonable  correlations.  In  the  initialization  and  prediction 
phase  of  the  present  study  the  heating  function  is  evaluated  according 
to  the  flow  chart,  figure  3,  which  is  self  explanatory. 

Large  scale  dynamical  prediction  of  tropical  disturbances  includ¬ 
ing  a  formal  parameterization  of  the  cumulus  convection  have  been  carried 
out  by  Manabe  and  Smagor insky  (1967).  It  would  be  pertinent  to  review 
their  parameterization  procedure.  In  their  general  circulation  model 
moist  convective  heating  and  condensation  takes  place  when  the  lapse 
rate  tends  to  become  super-moist  adiabatic  and  relative  humidity  tends 
to  exceed  100  per  cent.  During  this  situation  a  convective  adjustment 
is  forced  upon  the  numerical  solution  in  such  a  way  as  to  make  the  re¬ 
lative  humidity  and  the  equivalent  potential  temperature  uniform  in  the 
unstable  layer.  Furthermore  this  adjustment  is  carried  out  in  such  a 
way  as  to  keep  the  sum  of  potential,  internal,  and  latent  energy  invari¬ 
ant.  One  of  the  obvious  limitations  of  this  procedure  is  that  no  con¬ 
vective  process  on  small  scale  is  permitted  if  relative  humidity  <  100 
per  cent.  Synoptic  experience  in  the  tropics  is  quite  to  the  contrary 
(e.g.  Riehl,  19^5;  Baumhefner,  1968).  Manabe  and  Smagorinsky  (1967) 
discuss  a  rather  interesting  tropical  disturbance  that  forms  in  their 
long-term  general  circulation  experiment.  Its  large-scale  structure  is 
quite  similar  to  that  of  a  weak  tropical  storm,  although  a  large  inner 
core  of  this  disturbance  is  very  moist.  The  parameterization  procedure 
that  has  been  adopted  here  permits  convective  heat  and  moisture  flux  in 
the  vertical  direction  even  when  the  relative  humidity  is  less  than  100 
per  cent.  However,  this  scheme  presently  suffers  frem  the  lack  of  a 
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built  in  check  on  convective  adjustment  which  would  control  vertical 
soundings  and  force  them  to  look  like  their  atmospheric  counterparts. 

Such  a  built-in  convective  adjustment  is  considered  essential  for  long- 
term  general  circulations  experiments.  In  this  experiment  calculations 
were  extended  by  time  steps  to  U8  hours  and  no  severe  problems  such  as 
supersaturations  or  super-moist  adiabatic  instabilities  were  found. 

Kuo’s  scheme,  used  in  this  study,  behaves  roughly  in  the  following  manner 

A  conditionally  unstable  region  initially  experiencing  a  net 
convergence  of  moisture  flux,  I  >  0,  is  modified  in  a  matter  of  a  few 
hours  so  that  large-scale  values  of  T  approach  Tg,  and  q  approaches  qg. 
This  causes  "a"  to  become  smaller  and  the  region  to  become  sufficiently 
moist  that  further  mass  flux  into  the  region  begins  to  bring  in  relative¬ 
ly  drier  air  and  I  changes  sign  to  I  <  0.  When  this  happens  we  cease  to 
parameterize  the  convective  scale. 

(c)  Evaporation  and  precipitation 

Conservation  of  water  vapor  is  given  by  equation  (5).  Where  M 
represents  the  net  sources  and  sinks  of  water  per  unit  mass  of  air. 

Rohwer  (l93l)  expresses  evaporation  from  the  ocean  surface,  E  , 

v 

by  the  relation: 

Ev  =  (2.6  ♦  7.7  V10Q0)  (0.98  eu  -  e1000>/(2**-  *  3600 >.  (38) 

2 

The  units  of  are  gm/(cm  sec);  V^ooo  Is  1000  millibar  wind 
in  meters/sec;  e^  is  the  saturation  vapor  pressure  in  millibars  at  the 
temperature  of  the  water’s  surface;  and  e  is  the  vapor  pressure  of 
air  at  the  1000  millibar  surface.  No  attempt  was  made  to  vary  the  rough¬ 
ness  constant  for  variable  wind  speeds  over  the  oceans.  Over  land  areas 
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was  set  equal  to  zero. 


Webb  (i960)  has  critically  reviewed  the  use  of  such  an  empirical 
formula  for  estimating  evaporation  and  states  that  studies  over  Lake 
Eucumbe  show  that  error  would  be  within  10  per  cent. 

Precipitation  was  assumed  to  be  entirely  convective  and  was  evalu¬ 
ated  at  each  layer  using  Kuofs  parameterization  procedure,  figure  3. 


K 


C  a  (T  -  T) 
P _ S _ 

LAt 


(39) 


K  refers  to  the  vertical  level  and  P  has  the  units  sec 

K 

In  our  calculations  the  atmosphere  is  divided  into  four  layers, 
1000-800,  800-600,  600-U00  and  kOO  to  200  millibars.  Evaporation  is  only 
permitted  in  the  lowest  layer,  and  M  is  accordingly  defined  by  the 
relation : 


M 


E  C  a  (T  -  T) 

V  p  o _ 

g  Ap  ”  LAt 


(1*0) 


with  M  =  - 


C  a  (T_  -  T) 

E§ 


LAt 

of  the  lowest  layer. 


in  the  other  three  layers.  Ap  is  the  thickness 


5.  Surface  friction  and  diffusion  terms  of  model 

These  terms,  F^  and  F^,  in  the  momentum  equations  contain  two 
parts  in  this  model: 

i)  Subgrid  scale  lateral  and  vertical  diffusion  of  momentum  and 
ii)  effect  of  surface  wind  stress. 

This  formulation  is  almost  similar  to  that  used  by  Mintz  ( 196U )  in  his 
general  circulation  studies. 


We  write : 
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o  o  a2  3t 

n  ^  piC  O  U.  X 

Fx  =  un,  V  u  *  k  +  g  — 

3p 


and 


o  o  a2,,  9T,r 

Fy  =  7  v  +  <  “a  +  Zjf 

9p 


(Ul) 


(1*2) 


U  is  a  lateral  eddy  diffusion  coefficient  whose  value  in  the  present 

U  2-1 

experiment  is  5  x  10  meter  sec  and  k  is  a  vertical  eddy  diffusion 

-2  2  -I 

coefficient  whose  value  in  the  present  experiment  is  2  x  10  mb  sec 

The  primary  role  of  sub-grid  scale  lateral  and  vertical  diffusion  of 

momentum  is  to  remove  energy  from  the  smallest  permitted  scales.  It 

cannot  represent  physical  subgrid  scale  processes  in  the  horizontal  or 

vertical,  such  as  mesoscale  horizontal  jets  or  cumulus  convection.  It 

is  found  that  absence  of  vertical  diffusion  k  gives  rise  to  rather  large 

horizontal  wind  speeds  at  the  lowest  levels  within  a  ^8-hour  forecast 

period.  Calculations  did  not  appear  to  be  critically  dependent  on  the 

choice  of  the  lateral  coefficient. 

Surface  frictional  stresses  t  and  t  are  defined  in  terms  of 

x  y 

a  surface  roughness  constant  applied  at  the  lower  level,  p  =  1000  mb. 


2 

=  Cnpi/u'  +  v 


u 


(i*3) 


T 

y 


u2  + 


2 

v  v 


(kk) 


Cp  is  the  drag  coefficient,  whose  value  in  the  present  experiment  is 
2.5  X  10-3. 

Surface  stresses  prescribed  in  this  manner  provide  a  crude 
boundary  layer  for  our  model  between  the  1000  and  900  millibar  surfaces. 
Vertical  motions  at  the  900  millibars  are  influenced  significantly  by 
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the  presence  of  this  effect.  Frictional  convergence  of  mass  enhances 
low-level  convergence  of  moisture  and  the  latter  prescribes  a  percent 
area  of  convective  cloud  distribution  in  our  model.  A  discussion  of 
these  aspects  is  presented  with  the  results  of  the  numerical  experiment. 

6.  A  description  of  the  Initial  State:  March  1,  1965 >  1200  Z. 

It  is  not  possible  to  describe  all  of  the  scalar  fields  that 
define  the  initial  state.  A  selection  of  some  of  the  pertinent  fields 
is  presented  here. 

a)  Streamlines  of  the  observed  wind  field  at  1000  and  200  millibar 

surfaces .  Figures  4a  and  4b,  respectively,  portray  the  streamlines 
(solid  lines)  and  isotachs  (dashed  lines)  at  these  two  levels. 

The  following  synoptic  features  of  the  intertropical  conver¬ 
gence  zone  may  be  noted  at  the  lower  level.  An  asymptote  of  conver¬ 
gence  extends  from  152°W  to  about  l44°E  near  5°  north  latitude.  The 
speed  field  along  this  asymptote  is  quite  variable.  Wind  speed  ranges 
from  5  to  25  knots  over  most  of  the  map.  There  appears  to  be  no  simple 
structure  for  the  isotach  field  in  the  vicinity  of  this  asymptote  of 
convergence.  Wind  speed  maxima  can  be  found  along  and  north  of  this 
asymptote. 

There  is  a  belt  of  easterlies  over  most  of  the  northern  hemi¬ 
spheric  portion  of  the  analysis.  An  asymptote  of  convergence  near 
108°E  is  shown  as  a  part  of  a  cyclonic  eddy  at  3°N  near  l40°E.  In  the 
southern  hemispheric  tropics  interesting  features  are  the  two  cyclonic 
disturbances  at  1T°S,  110°E  and  at  15°S,  129°E.  The  speed  field  is 
shown  to  be  somewhat  organized  in  these  disturbances.  Over  the  rest  of 
the  Southern  Hemisphere  there  are  no  well  marked  centers. 
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At  200  mb  there  are  two  regions  of  westerly  wind  maxima,  one 
in  the  northwestern  portion  of  the  map  over  eastern  Asia,  and  another 
over  central  Australia  in  the  Southern  Hemisphere.  Over  the  equatorial 
latitudes  a  broad  easterly  current  with  speeds  near  1+0  knots  prevails . 

Four  anticyclonic  circulation  centers  (two  in  each  hemisphere)  consti¬ 
tute  an  important  part  of  the  upper-level  flow.  There  are  several 
other  transient  centers  that  are  portrayed  in  the  streamline  analysis 
of  this  chart. 

b)  Moisture  analysis  at  900  millibar  surface:  Moisture  analysis  was 
carried  out  by  first  using  all  available  stations  at  the  surface  level. 

A  satellite  nephanalysis  (valid  at  approximately  this  time)  was  used  as 
an  underlay  and  the  surface  streamline  features  were  used  to  complete 
the  subjective  distribution  of  relative  humidity  at  1000,  850,  700  and 
500  millibar  surfaces.  The  upper  levels  (850,  700  and  500)  were  arrived 
at  by  subjective  examination  of  cloud  cover,  upper  level  data  and  the 
surface  analysis.  The  relative  humidity  does  not  exceed  saturation 
values  on  the  large  scale  at  any  point.  The  analysis  is  very  subjective 
and  the  available  information  is  merely  a  guide.  The  distribution  of 
the  moisture  variable  is  very  poorly  defined,  hence  it  is  somewhat  unre¬ 
presentative.  The  initial  values  of  specific  humidity  at  900  (interpo¬ 
lated  between  1000,  850  millibars)  are  shown  in  figure  5(a).  A  moisture 
maximum  is  found  north  of  Australia  in  the  equatorial  Pacific  Ocean. 

It  has  values  of  specific  humidity  >  l8  gm/kgm.  Another  moisture  maxima 
center  with  values  close  to  20  gm/kgm  lies  northeast  of  Australia  and 
just  south  of  the  equator.  Over  the  continent  of  Australia  air  is 
fairly  dry,  values  being  close  to  5  gm/kgm.  Over  most  of  the  equatorial 
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tropics  values  are  near  15  gm/kgm.  There  is  a  gradual  decrease  of 
moisture  between  l6°N  and  28°N  to  around  6  gm/kgm  near  28°N. 

Because  of  the  uncertainty  in  the  analysis  of  the  moisture  vari¬ 
able  a  prediction  experiment  was  also  carried  out  with  uniform  initial 
relative  humidity  at  each  pressure  surface.  This  experiment  is  dis¬ 
cussed  separately  below. 

c)  Temperature  distribution  in  the  middle  troposphere:  500  millibar 
temperatures  are  shown  in  figure  5(b).  Thermal  structure  over  the 
central  Pacific  is  characterized  by  two  nearly  isothermal  warm  centers 
with  temperature  around  269°K.  These  warm  regions  define  the  warm  core 
structure  of  the  anticyclones  that  extends  to  the  200  millibar  surface; 
see  figure  Mb). 

In  the  northern  and  southern  portions  of  the  maps  two  well- 
marked  baroclinic  zones  are  present .  They  define  regions  where  wind 
maxima  are  found  at  the  200  millibar  surface,  figure  4(b).  The  active 
portion  of  the  intertropical  convergence  zone  lies  in  the  vicinity  of 
3°N  and  lU0°E.  In  the  middle  troposphere  mesoscale  temperature  anomalies 
are  likely  to  be  present  over  this  region.  The  smooth  large-scale  iso¬ 
thermal  patterns  could  be  interpreted  either  as  a  lack  of  smaller 
scale  data  or  an  efficient  lateral  and  vertical  mixing  process.  This 
could  be  possible  if  vertical  momentum  exchange  in  the  convective  re¬ 
gion  reduces  the  vertical  wind  shear  (Gray,  1967). 

Sea  surface  temperatures  were  obtained  from  the  monthly  mean 
charts  produced  by  U.S.  Bureau  of  Commercial  Fisheries  (1962)  and 
Hoyskoi  Atlas  (1953)  and  were  kept  fixed  during  the  experiment.  Ter¬ 
rain  height  distributions  were  obtained  from  the  tabulations  of  Berkofsky 
and  Bertoni  (i960). 
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d)  Satellite  cloud  cover:  During  this  period  (March  1-3,  1965)  TIROS  X 
made  a  few  passes  and  relevant  pictures  of  cloud  cover.  Figure  6  (a  and 
b)  shows  the  cloud  cover  from  approximately  130°E  to  l80°E  and  between 
20°N  and  20°S.  The  intertropical  convergence  zone  is  very  well  marked 
with  bright  clouds  near  North  between  l60°E  and  l8o°E  on  March  1. 

It  is  somewhat  more  extended  and  brighter  on  the  2nd  and  again  well 
marked  at  about  the  same  latitude  on  the  3rd.  In  the  vicinity  of  170°E 
and  15°S  a  well-marked  cloud  band  oriented  northwest-southeast  is  observ¬ 
able  on  all  three  days.  On  the  1st  and  2nd  a  bright  band  of  cloud  cover 
is  also  persistent  in  the  region  130°E  to  ll*0°E  and  10°S  to  20°S.  This 
is  over  north  central  and  north  western  Australia.  The  three  above-men¬ 
tioned  regions  of  intense  cloud  cover  are  found  in  the  vicinity  of 
cyclonic  surface  systems  and  asymptotes  of  low-level  convergence,  figure 
U(a).  Over  other  regions  the  amount  of  cloudiness  is  variable  during 
the  three-day  period.  The  picture  time  in  each  case  is  within  6  hours 
of  the  synoptic  time  of  analysis. 

No  attempt  will  be  made  to  relate  each  cloud  feature  observed  on 
these  TIROS  passes  with  results  of  computations.  Such  an  attempt  would 
be  futile  at  the  present  state  of  our  work.  Some  reference  will  be  made 
to  the  observed  cloud  cover,  at  least  in  its  larger-scale  context,  in 
the  following  sections. 

Me  shall  next  present  some  of  the  results  of  the  initialization 
from  the  balance  equations. 

e)  Initial  streamfunction  and  vorticity  fields:  Results  are  presented 
for  the  200  and  the  1000  millibars  surfaces,  figures  7(a  and  b),  respec¬ 
tively.  The  solid  isopleths  in  these  figures  show  the  streamfunction 
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and  the  dashed  isopleths  show  the  distribution  of  absolute  vorticity 
(units  10  ^  sec  ^). 

The  streamfunction  distribution  at  200  mb  contains  one  large 
anticyclonic  circulation  center  in  the  northern  Pacific  Ocean  and  a 
center  over  northern  Australia.  Strong  easterlies  are  found  over  the 
equatorial  latitudes  between  these  centers.  At  higher  latitudes  an  upper 

troughs  in  the  westerlies  contain  values  of  absolute  vorticity  ^1  x  10 

-1  -4  -1 

sec  in  the  Northern  Hemisphere  and  ^  -  1  x  10  sec  in  the  Southern 

Hemisphere.  The  zero  isopleth  of  absolute  vorticity  meanders  across  the 
equator.  At  the  1000  millibar  surface  a  maximum  value  of  absolute  vor¬ 
ticity  is  found  southeast  of  the  Truk  Islands  in  the  northern  Pacific 
Ocean.  This  is  associated  with  a  cyclonic  disturbance  along  the  western 
end  of  the  asymptote  of  convergence.  Observed  wind  directions  and  the 
streamfunction  are  nearly  parallel  at  the  200  millibar  surface,  while 
they  are  not  nearly  so  in  the  vicinity  of  the  asymptote  of  convergence 
at  the  1000  millibar  surface.  A  stronger  divergent  component  of  wind 
was  present  at  the  lower  level  in  the  vicinity  of  the  ITCZ.  The  cyclonic 
disturbances  over  the  subtropical  latitudes  of  the  Southern  Hemisphere 
are  indicated  by  the  large  negative  values  of  absolute  vorticity. 
f )  Vertical  Motions  on  top  of  the  friction  layer  and  calculated  cloud 

cover:  Figure  8  shows  the  isopleths  of  the  initial  balanced  vertical 
motions  x  10  ^  mb/sec  at  the  900  millibar  surface  (solid  lines),  while 
the  calculated  percent  area  of  convective  cloud  cover  is  shown  by  dashed 
lines.  The  shaded  region  has  some  fraction  of  the  synoptic  scale  grid 
(2°  latitude  x  2°  longitude)  covered  by  convective  clouds,  and  the  un¬ 
shaded  region  is  free  of  convective  clouds.  The  subtropical  disturb¬ 
ances  of  the  Southern  Hemisphere  exhibit  marked  convective  cloudiness. 
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Maximum  active  cloud  amounts  range  from  0.5  to  1  per  cent  in  these  dis¬ 
turbances.  Balanced  vertical  motion  in  these  disturbances  has  a  maximum 
value  os-  60  x  10~^  mb/sec  which  is  o^l  cm/sec.  The  principal  region 
of  interest  in  this  study  is  the  surface  asymptote  of  convergence.  As 
stated  earlier,  the  validity  of  the  balance  vertical  motions  near  the 
equator  are  questionable  because  Rossby  number  ^1.0. 

A  wide  region  of  maximum  rising  motion  (-  60  x  10  ^  mb/sec)  with 
maximum  cloud  amounts  ^1  per  cent  is  found  in  the  vicinity  of  the  ITCZ. 
Satellite  cloud  cover  shows  a  narrow  band  of  clouds,  figure  6* (a  and  b), 
along  U°N  near  170°E  and  a  region  with  widespread  convective  activity 
between  1^0°E  and  150°E.  A  region  of  sinking  motion  extends  from  136°E 
to  152°E  along  the  equator,  this  implies  that  air  reaching  the  equator 
from  the  Southern  Hemisphere  is  sinking  at  the  rate  of  1  cra/sec  before 
it  begins  to  ascend  in  the  vicinity  of  the  asymptote  of  convergence  in 
the  Northern  Hemisphere.  These  calculations  of  cloud  cover  and  the  in¬ 
tensity  of  vertical  motion  appear  somewhat  unrealistic  in  the  vicinity 
of  the  equator.  Short-range  prediction  with  this  initial  distribution 
of  mass,  motion  and  moisture  variables  reveals  some  interesting  changes 
in  the  vicinity  of  the  equator;  this  is  discussed  in  the  following 
section. 

7.  A  discussion  of  the  time  evolution  of  various  fields 

Some  of  the  results  of  the  calculations  will  be  discussed  in 
this  section. 

a)  1000  and  200  millibar  wind  field  at  2k  hours : 

Figures  9  (a  and  b),  respectively,  show  these  fields.  The  plotted 
winds  are  the  observations  at  2k  hours,  and  the  isopleths  are  the  predicted 
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are  extremely  interesting.  The  wind  field  at  1000  millibars  exhibits 
an  east-west  oriented  asymptote  of  convergence.  Convective  activity 
is  rather  well  marked  along  this  asymptote  of  convergence.  We  felt  that 
the  forecasts  after  30  hours  were  not  very  reliable,  mainly  because  of 
boundary  errors.  The  major  asymptote  of  convergence  and  associated  con¬ 
vective  cloud  band  exhibits  an  equatorward  movement.  Convective  activ¬ 
ity  weakens  somewhat  as  the  band  moves  from  U°N  at  2k  hours  to  near  the 
equator.  At  30  hours  this  system  strengthens  again  and  is  located  at 
5°S  at  U8  hours  of  forecast  time.  Meanwhile  a  new  asymptote  of  converg¬ 
ence  forms  at  12  to  l6°N  latitudes.  Our  calculations  show  that  rather 
large  changes  can  occur  in  the  tropics  during  a  ^8-hour  period.  This  is 
not  surprising  because  satellite  photographs,  figure  6  (a  and  b) ,  indi¬ 
cate  marked  changes  in  cloud  cover.  Manabe  and  Smagorinsky  (1967)  also 
mention  that  tropical  disturbances  have  a  rather  variable  structure 
during  1  to  3  day  periods. 

c )  Vertical  motions  at  900  millibars  (10  ^  millibars/sec) 

During  the  course  of  the  time  integration  other  relevant  fields 
like  convergence  of  mass  between  1000  and  900  millibars,  convergence  of 
moisture  between  1000  and  900  millibars  and  the  net  convergence  of 
moisture  between  1000  and  100  millibars  were  printed  out  at  6-hour  in¬ 
tervals.  It  is  true,  that  numerical  models  cannot  be  expected  to  reveal 
any  more  physics  than  what  is  already  stated  by  the  basic  equations  of 
our  models.  We  did  not  impose  any  very  direct  relationships  in  our 
basic  equations  between  the  cloud  cover  computations  and  the  various 
parameterization  features  listed  above.  We  observed  the  following  re¬ 


sults  of  interest: 
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Asymptotes  of  convergence,  percent  cloud  cover,  upward  vertical 
motions  on  top  of  friction  layer  and  net  moisture  convergence  all  re¬ 
vealed  rather  similar  patterns.  Figure  11  (a,b,c,d)  shows  the  friction 
layer  vertical  motions  in  6-hour  steps  with  superimposed  1000  millibar 
wind  directions.  Results  at  18  and  2b  hours  again  appeared  very  realistic. 
Largest  vertical  motions  are  of  the  order  -  300  x  10  ^  mb/sec,  i.e.  2s  U  cms 
rising  motion. 

The  evolution  of  the  vertical  motion  field  in  a  vertical  plane  is 
shown  in  figure  12.  The  vertical  plane  is  a  north-south  cross  section 
from  28°N  through  28°S  at  152°E  longitude.  This  plane  intersects  an  act¬ 
ive  portion  of  the  ITCZ.  The  field  of  vertical  motion  undergoes  consider¬ 
able  adjustment  in  the  first  18  hours  in  the  equatorial  latitudes.  This 
field  appears  very  realistic  at  18  and  2b  hours  of  forecast  time.  The 
largest  magnitudes  of  vertical  motion  occur  at  2b  hours  .and  are  -  200  x 
10  mb/sec  (*  3  cm/sec)  and  +  500  x  10  mb/sec  (**  -  7  cm/sec)  in  the 
mid  troposphere.  The  adjustment  of  the  vertical  motion  field  produces  a 
solution  that  is  quite  asimilar  to  initial  vertical  motions.  Magnitudes 
of  “10qq  are  of  the  order  of  25  x  10-^  mb/sec  (vertical  motion  is  zero  at 
100  millibars). 

d )  Some  parameters  related  to  parameterization  of  cumulus  convection : 

In  figure  13  we  show  the  tine  variation  in  6-hour  steps  at  152°E 
and  between  28°S  and  28°N  of  the  following  quantities : 

i)  Specific  humidity  at  900  mb  (gm/kgm) 

ii)  Percent  cloud  cover 

8  1 

iii)  Net  moisture  convergence  (unit  10  mb  meter  sec). 
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The  field  of  specific  humidity  changes  very  slightly  during  the 
first  2k  hours.  Along  this  section  specific  humidity  varies  between  10 
and  18  gm/kgm  between  the  equatorial  maximum  value  and  lower  values  near 
28°N  and  28°S.  The  region  of  maximum  moisture  is  not  the  most  convec- 
tively  active  region.  Percent  cloud  cover  (maximum  1  per  cent)  and  the 
net  moisture  convergence  are  very  nearly  proportional  to  each  other. 

These  fields  undergo  rapid  variations  during  the  first  18  hours. 

It  appears  from  the  nature  of  these  variations  that  spurious  ver¬ 
tical  advections  of  heat  and  moisture  take  place  during  the  adjustment 
process.  For  short-range  predictions  in  low  latitudes,  it  is  highly 
desirable  that  such  spurious  effects  be  removed  right  from  the  start. 
Unfortunately ,  the  present  model  is  incapable  of  alleviating  this  prob¬ 
lem.  To  the  best  of  the  author's  knowledge,  suppression  of  this  initial 
oscillation  has  not  been  satisfactorily  carried  out  in  any  real  initial 
data  experiments.  Miyakoda  and  Moyer  (1968)  and  Nitta  and  Hovermale 
(1967)  have  proposed  procedures  for  initialization  that  could  alleviate 
this  difficulty,  but  because  of  limited  tests  in  their  studies,  no 
general  statement  can  be  made  at  present. 

e)  An  experiment  with  constant  relative  humidity  on  each  pressure  level 

at  initial  time : 

As  stated  earlier,  the  initial  moisture  distribution  is  somewhat 
unreliable,  hence  relative  humidity  obtained  from  subjective  analysis  was 
averaged  horizontally  over  the  domain  of  integration  in  this  experiment. 
Thus  initially  no  gradients  of  relative  humidity  were  present  on  constant 
pressure  surfaces  but  it  should  be  noted  that  specific  humidity  is  not  a 
constant.  The  results  of  numerical  prediction  of  the  mass  and  motion 
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fields  at  2U  hours  were  quite  similar  to  those  obtained  from  detailed 
analysis  of  moisture.  This  implies  the  following  important  feature  of 
the  parameterization  of  the  cumulus  scale  in  this  experiment: 

V  3  q\V  qV  »  W,  i.e.W  3  V  q 

is  small.  Mass  convergence  is  a  more  pronounced  feature  of  the  flows 
and  accounts  for  the  net  convergence  of  flux  of  moisture  while  hori¬ 
zontal  advection  of  moisture  contributes  less.  This  statement  must  be 
interpreted  with  seme  caution,  because  this  is  only  true  in  a  gross 
sense.  The  author  has  carried  out  another  detailed  synoptic  study 
(fCrishnamurti ,  1968b)  of  various  parameters  of  large-scale  flows  and 
their  relation  to  convective  cloud  cover  as  revealed  by  radar  returns 
over  the  Caribbean.  It  was  found  that  mass  convergence  and  moisture 
convergence  were  not  of  the  same  sign  over  all  regions  near  900  mb. 

Since  the  moisture  variable  is  one  of  the  most  unreliable  fields  in 
our  input  we  might  suggest  that  a  constant  distribution  of  relative 
humidity  on  pressure  surfaces  might  be  a  useful  compromise  for  very 
short-range  prediction. 

8.  Some  comments  on  the  dynamics  of  the  intertropical  convergence  zone 
The  field  of  vertical  motions  in  the  balance  initial  state  is 
in  good  agreement  with  the  synoptic  disturbances  north  of  15°N  and 
south  of  15°S.  In  the  vicinity  of  the  equator  initial  vertical  motion 
field  shows  a  general  region  of  rising  motion  north  of  the  surface  asymp¬ 

tote  of  convergence  and  does  not  appear  to  be  very  realistic.  A  short 
range  prediction  in  a  matter  of  18  to  2h  hours  yields  a  rather  realistic 


pattern  of  rising  motion  that  extends  from  east  to  west  almost  along 
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the  asymptote  of  convergence  in  the  surface  forecast  wind  field.  The 
initial  discrepancy  must  be  due  to  the  inadequate  dynamics  of  the  bal¬ 
anced  initial  state  which  is  only  valid  for  flows  characterized  by  a 
small  Rossby  number.  The  adjustment  of  the  three-dimensional  motion 
field  to  the  temperature,  moisture,  mass  field,  and  other  external 
boundary  parameters  is  interesting.  A  narrow  axis  of  computed  cloud 
cover  at  2b  hours  along  ^°N  raises  the  question  whether  the  output  in¬ 
formation  at  2h  hours  could  be  used  to  investigate  some  aspects  of  the 
dynamics  of  the  intertropicaJ  convergence  zone.  A  short-range  predic¬ 
tion  is  certainly  not  a  simulation  experiment,  our  initial  data  does 
already  contain  the  proper  ingredients  in  the  motion,  temperature  and 
moisture  variables  to  resolve  in  a  matter  of  2b  hours  a  narrow  region 
of  cloud  cover  and  rising  motion  along  the  asymptote  of  low-level  con¬ 
vergence.  At  best,  we  could  hope  to  merely  carry  out  a  diagnostic  study 
with  the  detailed  output  information  on  tape  at  2b  hours.  At  this  time 
there  is  some  sort  of  a  balance  of  the  dependent  variables  for  large 
Rossby  number  flows  in  the  vicinity  of  the  equator. 

Diagnostic  studies  are  generally  limited  in  scope,  if  calcula¬ 
tions  are  carried  out  in  a  small  domain  (as  they  are  in  our  studies), 
any  description  of  the  local  circulation  can  be  somewhat  biased.  A  per¬ 
tinent  question  that  could  be  raised  at  this  point  would  be  on  the 
maintenance  of  the  kinetic  energy  of  the  wind  field. 

Let  us  consider  a  box  enclosing  the  region,  l6°N  through  l6°S, 
1000  millibars  through  100  millibars  and  from  10U°E  through  128  °U. 
Because  of  the  imposed  cyclic  continuity  there  are  no  boundary  fluxes 


in  the  zonal  direction. 
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The  following  integrals,  respectively,  represent  the  conversions 
of  zonal  kinetic  energy  and  eddy  available  potential  energy  into  eddy 
kinetic  energy,  see  Lorenz  (i960),  and  Manabe  and  Smagorinsky  (1967) • 


Equations  (**5)  and  (1*6)  represent  contributions  to  energy  transforma¬ 
tions  at  a  fixed  latitude,  primes  are  departures  from  zonal  averages 
and  bar  is  a  zonal  average,  other  symbols  are  explained  elsewhere. 

The  net  energy  conversion  over  the  box  may,  if  desired,  be  obtained  by 
merely  integrating  these  conversions  over  the  surface  area  of  the  box. 
Both  of  these  quantities  may  be  expressed  in  millibar  meter  sec  ^  units. 
This  unit  seems  more  appropriate,  because  we  are  merely  comparing  the 

two  rates  of  energy  transfer  with  respect  to  each  other.  For  similar 

2 

calculations  Manage  and  Smagorinsky  (1967)  use  Joule/(cm  day)  as  the 

2 

unit,  the  conversion  between  these  sets  of  units  is  1  Joule/(cm  day)  = 
(g^O  mb  meter  sec-1. 

Figure  lh  shows  the  results  of  calculations  at  t  =  2**  hrs.  In 
the  vicinity  of  the  intertropical  convergence  zone,  the  so-called  baro- 
tropic  conversion  is  a  transfer  from  eddy  to  zonal  kinetic  energy  while 
there  is  a  stronger  conversion  of  eddy  potential  to  eddy  kinetic  energy. 
This  result  is  not  surprising  if  we  examine  the  synoptic  maps  at  differ¬ 
ent  levels  at  t  =  2U  hours. 

The  wind  field  at  1000  and  200  millibars  exhibits  a  belt  of 
easterlies  near  the  ITCZ.  Easterlies  are  strong  at  low  levels  and  at 
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higher  levels,  being  ouite  variable  and  weak  at  the  600  and  1*00  millibar 
surfaces.  At  the  lower  level  intense  rising  centers  u'  <0  are  found 
where  u'  >  0  (weaker  easterlies  near  the  ITCZ,  see  figure  15)  in  this 
region  ^  <0.  At  higher  levels  w'  <  0,  u'  <0  and  >  0.  The  net  re- 
suit  of  these  configurations  of  the  motion  field  yields 
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Furthermore ,  along  the  asymptote  of  convergence  lies  a  region  of  weak 

surface  winds,  figure  15.  Just  north  of  the  ITCZ,  along  a  latitude 

circle  v'  <  0,  u'  >0,  ^  <  0,  south  of  the  ITCZ  v'  >  0,  u'  >  0,  >  0, 

a,v  ay 

and  hence 
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Thus  ,  the  geometric  configuration  of  the  computed  velocity  field  at  2h 

hours  has  a  net  barotropic  energy  exchange  (K^  •  K^,)  <  0.  This  exchange 

is  more  than  compensated  by  the  ascending  warm  air  (P^  °  K^)  near  the 

ITCZ,  figure  11.  The  barotropic  exchange  away  from  the  ITCZ  shows 

(K^  •  Kg)  >  0  in  both  the  Northern  and  Southern  Hemispheres,  this  result 

is  in  accordance  with  Charney's  (1963)  scale  analysis  for  large-scale 

motions  in  the  tropics.  The  dominant  energy  conversion  term  is  (P^  •  K^,), 

though  between  8°S  and  l6°S  and  between  8°N  and  l6°N  (K  •  K  )  is  an 

L  h 

appreciable  fraction  of  the  total  exchange. 

The  intense  rising  motion  (u)’)  along  the  ITCZ  in  our  model  is 
mainly  attributable  to  the  coupling  of  strong  frictional  convergence  and 
the  convective  release  of  latent  heat.  An  experiment  was  also  performed 
without  latent  heat  but  including  surface  friction.  The  results  indicated 
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a  breakdown  of  the  elongated  asymptote  of  convergence  into  several  large 
eddies  in  a  matter  of  18  hours.  Perhaps  this  result  was  a  consequence 
of  improper  energetics  and  the  adjustment  of  our  balance  initial  state 
to  this  situation. 

Because  of  the  limited  nature  of  the  present  investigation,  a 
complete  energetics  of  an  open  system  containing  the  ITCZ  was  not  carried 
out. 


10.  Some  concluding  remarks: 

The  procedure  for  parameterization  of  cumulus  convection  used  in 
this  study  is  rather  crude.  Vertical  momentum  transfer  in  sub-grid  scale 
motion  is  not  included.  Recently  Ooyama  ( 19 67)  and  Estoque  (1967)  have 
considered  this  effect  in  their  studies.  Its  inclusion  would  be  a  simple 
extension  of  the  present  study.  We  have  shown  that  a  somewhat  realistic 
short-range  prediction  of  larger  scale  flows  in  equatorial  latitudes  is 
not  possible  without  a  parameter is at ion  of  the  cumulus  scale.  We  have, 
however,  made  a  tacit  assumption  that  cumulus  convection  is  always  express¬ 
ible  as  a  function  of  larger-scale  moisture  convergence.  This  is  generally 
the  case,  but  it  excludes  the  possibility  of  isolated  cumulus  convection 
with  local  smaller-scale  mass  compensation.  Convection  over  heated  islands 
and  sea  breeze  effects  generally  occur  without  bearing  any  relationship  to 
organized  larger-scale  mass  or  moisture  convergence,  though  they  may  be 
modified  by  the  prevailing  synoptic  situation.  This  is  not,  however,  a 
hopeless  problem.  Parameterization  of  cumulus  convection  could,  in  prin¬ 
ciple,  be  carried  out  to  include  both  of  these  effects  separately.  There 
is  thus  an  obvious  need  for  a  detailed  climatological  investigation  of 
cloudiness  over  all  islands  and  coastal  regions  that  are  known  to  exhibit 
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local  diurnal  convection.  If  such  statistics  were  accumulated  over  the 
entire  tropics  then  this  information  could  be  introduced  in  some  empiri¬ 
cal  manner  into  larger  scale  numerical  prediction  models.  This  is  neces¬ 
sary  because  convective  clouds,  no  matter  how  they  are  formed,  do  transfer 
large  amounts  of  moisture,  heat  and  momentum  in  the  vertical  direction  and 
the  larger-scale  motions  are  gradually  altered  as  a  result.  If  a  coarse 
mesh  of  the  order  of  few  hundred  kilometers  is  used  in  prediction  experi¬ 
ments  then  it  may  be  important  to  include  this  second  class  of  convective 
element.  Parameterization  of  cumulus  convection  requires  a  knowledge  of 
the  distribution  of  moist  unstable  regions  in  the  horizontal  and  in  the 
vertical.  A  knowledge  of  the  distribution  of  the  moisture  variable  is 
thus  quite  important,  yet  we  find  that  the  initial  analysis  of  this  field 
is  very  difficult  especially  from  the  currently  available  data  networks. 
Several  massive  research  efforts  in  the  tropics  are  currently  being  con¬ 
templated  in  the  United  States  and  other  countries.  It  would  be  highly 
desirable  to  examine  ways  to  obtain  detailed  three-dimensional  distribu¬ 
tion  of  moisture  variables  in  low  latitudes. 

Some  of  the  important  results  of  the  present  study  are  as  follows: 

The  dynamical  structure  of  the  tropical  atmosphere  in  equatorial 
latitudes  is  poorly  described  by  the  balance  model,  but  is  reasonably  well 
represented  if  a  short-term  numerical  prediction  is  carried  out. 

At  2h  hours  of  forecast  time  the  following  features  are  of  inter¬ 
est.  A  well  marked  surface  asymptote  of  convergence,  the  ITCZ,  is  found 
at  h°N.  Along  this  asymptote  an  axis  of  frictional  mass  and  moisture  con¬ 
vergence  is  found  with  a  maximum  computed  convective  cloud  cover.  This 
study  reveals  a  strong  coupling  between  these  three  features.  Along  the 
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equator,  sinking  motion  and  absence  of  clouds  is  indicated  by  the  compu¬ 
tations.  Furthermore  we  find  that  in  the  vicinity  of  the  ITCZ,  energy 
transformations  show  that  (P^,  •  K^,)  is  positive,  in  agreement  with  the 
results  of  Fanabe  and  Smagorinsky  (1967)  but  we  also  find  a  surprising 
result  that  (Kg  •  K^,)  is  positive.  We  have  shown  that  the  latter  result 
is  possible  because  of  the  detailed  structure  of  wind  field  in  low  lati¬ 
tudes.  Away  from  the  ITCZ  we  find  that  (K„  »  K  )  is  positive  implying  a 
barotropic  instability  of  the  zonal  current.  This  study  has  demonstrated 
a  limited  success  in  short-range  numerical  prediction  from  real  initial 
data  in  low  latitudes. 

Besides  the  difficulties  of  the  parameterization  of  the  cumulus 
scale  other  major  defects  of  this  study  are  related  to  lack  of  observa¬ 
tions,  contamination  of  forecasts  due  to  lateral  boundary  conditions  and 
the  presence  of  gravity  inertio-oscillations  during  the  ^initial  phase  of 
the  forecast  period.  Further  work  is  continuing  to  examine  in  detail 
some  of  these  problems. 
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Figure  Captions 


Figure  1.  Sketch  of  nine  points  on  a  pressure  surface  used  in  quasi- 
Langrangian  advection.  P  is  the  origin  of  a  parcel  at  time 
t  -  At,  and  Q  is  a  grid  point  where  parcel  arrives  at  time  t. 

Figure  2.  Vertical  staggeriing  of  variables  in  the  initialization  and 
prediction  phase. 


Figure  3.  Procedure  used  for  the  parameterization  of  cumulus  convection. 


Figure  U.  (a)  (top)  Streamlines  (solid  lines)  and  isotachs  (dashed 

lines ,  knots )  for  March  1 ,  1965 ,  1200Z ,  at 
1000  millibars,  t  =  0  in  our  experiment. 

(b)  (bottom)  Streamlines  (solid  lines)  and  isotachs  (dashed 
lines,  knots)  for  March  1,  1965,  1200Z,  at  200 
millibars,  t  =  0  in  our  experiment. 


Figure  5*  (a)  (top)  Specific  humidity  distribution  at  900  millibars, 

units  gm/kgm,  for  March  1,  1965,  1200Z,  t  =  0  in 
our  experiment. 

(b)  (bottom)  Temperature  distribution  at  500  millibars, 

units  °A,  for  March  1,  1965,  1200Z,  t  =  0  in  our 
experiment. 


Figure  6.  (a)  (top)  Satellite  mosaics  for  March  1,  2  and  3,  1965,  in 

the  vicinity  of  lUo°E  longitude. 

(b)  (bottom)  Satellite  mosaics  for  March  1,  2  and  3,  1965,  in 
the  vicinity  of  1T0°E  longitude. 
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Figure  7*  (a)  (top)  Initial  streamfunction  at  200  millibar  surface, 

units  10^  meter^  sec  ^ ,  March  1,  1965,  1200Z. 

(b)  (bottom)  Initial  stream  function  at  1000-millibar  surface, 
units  10^  meter^  sec  March  1,  1965,  1200Z. 

Figure  8.  Initial  balance  vertical  motions  (solid  lines)  at  900  millibar 
units  10^  millibar/sec.  The  initial  calculated  percent  convec¬ 
tive  cloud  cover  is  shown  by  dashed  lines  and  this  region  is 
shaded.  March  1,  1965,  1200Z. 

Figure  9-  (a)  (top)  24-hour  forecast  of  the  wind  field  at  1000  milli¬ 

bars,  solid  lines  are  streamlines  and  dashed 
lines  are  isotachs  in  knots.  The  plotted  winds 
are  observations  at  verification  time,  i.e.  March 
2,  1965,  1200Z. 

(b)  (bottom)  Same  as  figure  9(a)  except  at  200  millibars. 

Figure  10.  (a,  b,  c  and  d)  6-,  12- ,  18-  and  24-hour  forecast  of  the  dis¬ 
tribution  of  percent  area  of  convective  cloud 
cover,  shaded  region.  The  wind  arrows  are  the 
directions  of  the  forecast  wind  vectors  at  the 
1000  mb  surface. 

Figure  11.  (a,  b,  c  and  d)  6-,  12- ,  18-  and  24-hour  forecast  of  the  dis¬ 
tributions  of  900-millibar  vertical  motions, 
units  10~5  millibar/sec,  the  shaded  region  indi¬ 
cates  rising  motions.  The  wind  arrows  are  the 
directions  of  the  forecast  wind  vectors  at  the 


1000  mb  surface. 


Fig.  12. 


Figure  13 


Figure  l1* 


Figure  15 


Vertical  cross-section  of  vertical  motion  (10  ^  millibar/sec) 
at  6,  12,  18  and  2k  hours  of  forecast  time.  This  cross 
section  is  taken  in  a  south-north  plane  at  152°W  longitude, 
and  it  intersects  the  ITCZ. 

Time  variation  (along  ordinate)  of  selected  parameters  at 
152°W  from  south  to  north  (along  abscissa). 

q  specific  humidity  at  900  mb  gm/kgm 
a  percent  cloud  cover 

8  l 

I  net  moisture  convergence  (10  mb  meter  sec). 

Selected  energy  conversions  at  t  =  2k  hours  of  forecast  time, 
evaluated  from  equations  (1*5)  arid  (1*6)  of  text. 

An  enlarged  view  of  the  1000  mb,  2U-hour  forecast  of  wind  direc- 
tions( arrows)  and  wind  speeds  (knots,  solid  lines)  in  the 
vicinity  of  the  ITCZ  (heavy  dark  line). 
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Figure  1.  Sketch  of  nine  points  on  a  pressure  surface  used  in  quasi- 
Langrangian  advection.  P  is  the  origin  of  a  parcel  at  time 
t  -  At ,  and  Q  is  a  grid  point  where  parcel  arrives  at  time  t. 
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Figure  2.  Vertical  staggeriing  of  variables  in  the  initialization  and 


prediction  phase. 


FLOW  CHART  FOR  CUMULUS  SCALE  HEATING:- 


Figure  3*  Procedure  us^d  for  the  parameterization  of  cumulus  convection. 


Figure  4.  (a)  (top)  Streamlines  (solid  lines)  and  isotachs  (dashed 

lines,  knots)  for  March  1,  1965,  1200Z,  at 
1000  millibars,  t  =  0  in  our  experiment. 


(b)  (bottom)  Streamlines  (solid  lines)  and  isotachs  (dashed 
linos,  knots)  for  March  1,  1965,  1200Z,  at  200 


millibars ,  t  =  0  in  our  experiment 


our  experiment. 


units  °A,  for  March  1,  1965,  1200Z,  t  =  0  in  our 


experiment. 


Figure  6.  (a)  Satellite  mosaics  for  March  1,  2,  and  3,  1965  in  vicinity  of  140°E  longitude. 
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Figure  6.  (b)  Satellite  mosaics  for  March  1,  2,  and  3,  1965  in  vicinity  of  170°E  longitude. 
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Figure  7.  (a)  (top) 


Initial  streamfunction  at  200  millibar  surface, 
units  10^  meter^  sec  \  March  1,  1965,  1200Z. 


(b)  (bottom)  Initial  stream  function  at  1000-millibar  surface, 

5  2-1 

units  10  meter  sec  , 


March  1,  1965,  1200Z 
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Figure  8. 


Initial  balance  vertical 
units  1C5  millibar/sec. 


motions  (solid  lines)  at  900  millibar 
The  initial  calculated  percent  convec¬ 


tive  cloud  cover  is  shown  by  dashed  lines  and  this  region  is 


shaded.  March  1,  1965,  1200Z. 


28*N 


(top)  24-hour  forecast  of  the  wind  field  at  1000  millibars, 
solid  lines  are  streamliner  and  dashed  lines  are 
isotachs  in  knots*  The  plotted  winds  are  observations 
at  verification  time,  i.e.  March  2,  196 5,  1200Z. 
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(b)  (bottom)  Same  as  figure  9(a)  except  at  200  millibars. 


shaded  region*  The  wind  arrows  are  the  directions  of  the 
forecast  wind  vectors  at  the  1000  mb  surface. 
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Figure  li. 


(a,  b,  c,  &  d)  6-,  i2T,  18-;,  ana  24-hour  forecast  ofj.the  dis¬ 
tribution  of  900-millibar  vertical  motions ,  units  lO’“'>millibar/ 
sec,  the  shaded  region  indicates  rising  motions*  The  wind 
arrows  are  the  directions  of  the  forecast  wind  vectors  at  the 
1000  mb  surface. 
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Fig.  12.  Vertical  cross-section  of  vertical  motion  {10-^  millibar/sec) 
at  6,  12,  18  and  2U  hours  of  forecast  time.  This  cross 
section  is  taken  in  a  south-north  plane  at  152°W  longitude, 
and  it  intersects  the  ITCZ. 
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Figure  13*  Tim*  variation  (along  ordinate)  of  selected  parameter*  »t 
152  V  froa  south  to  north  (along  abscissa). 
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Figure  lU.  Selected  energy  conversions  at  t  =  2k  hours  of  forecast  time, 
evaluated  from  equations  (45)  add  (46)  of  text. 


Figure  15.  An  enlarged  view  of  the  1000  mb,  24-hour  forecast  of  wind  direc- 


tions(arrows)  and  wind  speeds  (knots,  solid  lines)  in  the 
vicinity  of  the  ITCZ  (heavy  dark  line). 
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13.  ABSTR  AC  T 


In  order  to  obtain  the  dynamical  structure  of  the  tropical  atmosphere  in 
equatorial  latitudes  a  short  range  prediction  experiment  is  proposed.  The  initial 
field  contains  the  inter tropica  1  convergence  zone  and  associated  disturbances  over 
western  Pacific  Ocean  during  March  1965.  A  complete  initial  state  is  constructed 
using  a  consistent  balance  system  of  equations.  The  procedure  involves  construc¬ 
tion  of  pressure,  temperature  and  vertical  motion  distributions  starting  from  an 
observed  rotational  part  of  the  wind  field.  It  is  shown  that  this  procedure, 
whose  validity  assumes  a  small  Rossby  number,  does  not  yield  a  realistic  structure 
of  the  vertical  motion.  A  short  range  prediction  with  a  multi-level  prediction 
model  yields  some  interesting  solutions  in  the  vicinity  of  the  inter tropica  1  con¬ 
vergence  zone.  During  the  first  18  hours  of  prediction  an  adjustment  of  the  motion 
and  the  mass  field  ensues  with  gravity  inertia  oscillations.  A  detailed  discussion 
of  some  dynamical  aspects  at  24  hours  is  presented.  An  important  feature  of  the 
model  is  a  parameterization  of  cumulus  convection  as  a  function  of  large-scale 
moisture  convergence.  The  role  of  cumulus  scale  heating  in  the  vicinity  of  the 
intertropica 1  convergence  zone  is  explored  by  carrying  out  experiments  with  and 
without  diabatic  heating. 
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